Hopping induced continuous diffusive dynamics below the non-ergodic transition 
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In low temperature supercooled liquid, below the ideal mode coupling theory transition tempera- 
ture, hopping and continuous diffusion are seen to coexist. We present a theory which incorporates 
interaction between the two processes and shows that hopping can induce continuous diffusion in the 
otherwise frozen liquid. Several universal features arise from nonlinear interactions between the con- 
tinuous diffusive dynamics (described here by the mode coupling theory (MCT)) and the activated 
hopping (described here by the random first order transition theory ). We apply the theory to real 
systems (Salol) to show that the theory correctly predicts the temperature dependence of the non- 
exponential stretching parameter, /3, and the primary a relaxation timescale, r. The study explains 
why, even below the ergodic to non-ergodic transition, the dynamics is well described by MCT. The 
non-linear coupling between the two dynamical processes modifies the relaxation behavior of the 
structural relaxation from what would be predicted by a theory with a complete static Gaussian 
barrier distribution in a manner that may be described as a facilitation effect. Furthermore, the 
theory explains the observed variation of the stretching exponent j3 with the fragility parameter, D. 



I. INTRODUCTION 

The glass transition is characterized by a number of 
interesting kinetic phenomena. Slow non-exponential re- 
laxation of time correlation functions over large time win- 
dows is one such important phenomenon. This relax- 
ation is usually well approximated by the stretched ex- 
ponential, Kohhausch- William- Watts (KWW) formula 
1], (f){t) = exp(-{t/T)f, with both (1 and r exhibit- 
ing non trivial temperature dependence. The origin of 
the stretching is usually attributed to the presence of 
dynamic heterogeneity in the system 0, yj. The tem- 
perature dependence of the relaxation time can be de- 
scribed by the the Vogel-Fulcher-Tamman (VFT) expres- 
sion, r = tvftg DT ° ' ( T-T ° ), where tvft is the high tem- 
perature relaxation time, T is the VFT temperature and 
D is the fragility index. The fragility index, D determines 
the degree of deviation from the Arrhcnius law. Exper- 
imental and theoretical model studies have shown that 
(3 and D are correlated [J, Q ■ The temperature depen- 
dence of r is also described by phenomenological MCT 
expression, r - [(T- T/^/T/**] -7 , (7 > 0). T[ lt is 
referred to as the MCT transition temperature. Above 
T/ % \ MCT is found to explain many experimental results 
[I, 0, B S [3, and below T/ 1 *, the MCT description of 
continuous diffusion breaks down. It is conjectured that 
this breakdown is due to the ergodic to non-ergodic tran- 
sition in the dynamics and below T/ lt activated dynamics 
becomes a dominant mode of transport. However in an 
elegant work, Burmer and Reichman 11] (BR) have re- 
cently shown that the idealized MCT breaks down at a 



much higher temperature, Tj?, which corresponds to the 
temperature, Tt,, where the landscape properties change. 
Kob et al. have shown that the structural MCT predicts 
the proper dynamics till very close to T/ 1 ' but needs 
static inputs to be calculated at a higher effective tem- 
perature [l4|. Wolynes and co-workers have shown that 
below the ergodic to non-ergodic transition but before 
the predicted crossover to the activated d yna mics, string 
or fractal excitations drive the dynamics [131 ]. 

Computer simulation studies predict the coexistence 
of continuous diffusion and hopping as mechanism of 
mass transport at temperatures much above T/ 1 ' [lj, [l5| . 
These studies show that individual hopping events are 
often followed by enhanced continuous diffusion or more 
hops of the surrounding atoms or molecules [lj,ll6|]. Sim- 
ulations have also shown that a single hopping event re- 
laxes the local stress [lj], hence it is expected that hop- 
ping events are followed by continuous diffusive dynam- 
ics. As mentioned earlier, from BR study we know the 
temperature where the ergodic to non-ergodic transition 
takes place (T® = Tt,) [ll|. However, there is no clear 
theoretical understanding of the dynamics in the range, 
T[ %t < T < T®. Furthermore, we need to understand, 
why below 7* the MCT still seems to explain the dy- 
namics rather well and finally, what happens at T/ 1 ' that 
leads to the breakdown of the MCT. 

In an earlier work it was shown that indeed below 
the ergodic to non-ergodic transition, both activated and 
continuous diffusion mechanisms are present. These two 
are nonlinearly coupled and may drive each other 17| ■ It 
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was shown that the full dynamics is a synergetic effect of 
continuous diffusive motion and activated dynamics [171 ] . 
Here we show that the above coupling can lead to quite 
interesting phenomena. First for the Salol system we find 
that, the theory predicts the correct temperature depen- 
dence of relaxation time. Following BR analysis [11 1 we 
consider that the ergodic to non-ergodic transition in the 



ideal MCT occurs at Tj? = 278/c. However, when we fit 
the relaxation timescale to the phenomenological MCT 
expression, the timescale shows an apparent divergence 
at a lower temperature, thus predicting T[ %t = 258K 
which is similar to that found in experimental studies 6]. 
Thus the study demonstrates that well below the ideal- 
ized MCT ergodic to non-ergodic transition temperature 
where non-mean field activated events are present, the 
dynamics can be well described by MCT. We find that, 
across T[ %t , although both the transport mechanisms are 
operative and there is no discontinuity in their motion, 
there is an abrupt or sudden rise in the contribution from 
the activated events. This leads to the apparent break- 
down of the MCT. 

Second, the theory predicts the correct temperature 
dependence of the stretching parameter j3. We find that 
the bare barrier height distribution that describes iso- 
lated hopping, gets renormalized due to the coupling of 
the hopping dynamics to the continuous diffusive motion. 
The renormalized distribution of the barriers, referred to 
as the dynamic barrier height distribution, is different 
from the initial assumption of distribution referred to as 
the static barrier height distribution. There is a prefer- 
ence for low barrier hopping events leading to facilita- 
tion effect. Third, due to facilitation effect, the theory 
predicts a correct relationship between the stretching pa- 
rameter, /3, and the fragility index, D. 

The organization of the rest of the paper is as follows. 
In the next section we describe the theoretical scheme. 
Section 3 contains results and discussions. Section 4 con- 
cludes with a discussion on the results. 



II. THEORETICAL FORMULATION 

The present work uses a scheme of calculation similar 
to the one presented earlier [lj} with a few modifications. 
For describing the diffusive motion we use P12 model of 
the MCT [18] and we also describe the activated hopping 
dynamics with a static barrier height distribution [5[. In 
describing the dynamics the coupling between the differ- 
ent wave- vectors are neglected and the contribution from 
the static and dynamic quantities are calculated at a sin- 
gle wavenumber q = q m (q m is the wave number where 
the peak of the structure factor appears) which is known 
to provide the dominant contribution. The MCT part 
of the intermediate scattering function, 4>MCT{t) is self 
consistently calculated with the full scattering function, 
4>(t). In writing the equation of motion we consider that 
there are distribution of hopping barriers in the system 
arising from the entropy fluctuation [5|. We know that 
in the absence of hopping motion the MCT dynamics be- 
low certain temperature is completely frozen, leading to 
a separation of time scale between the hopping dynamics 
and MCT dynamics. Thus for the frozen MCT dynamics, 
it appears that multiple relaxation channels (present due 
to different hopping barriers) open up and all these chan- 
nels act in parallel to relax the system. Thus when the 



contribution from these multiple hoppings are included in 
the total relaxation, the equation of motion of MCT part 
of the intermediate structure function, 4>MCT{t), can be 
written as, 

<pMcr{t) + 7</>mct(£) + Stf)4>MCT{t) 

+ Aing / dt'cf> MCT (t')y2<t> s hop (t')4>Mc T (t-t>) 

Jo 

+ X 2 n 2 Q I dt'[(l)MCT(t')y2<t> S ho P (t')} 2 4>MCT(t-t') 

Jo 
= (1) 

Eqfl] is essentially of the same general form as derived 
by Sjo gren and Gotze for the extended mode coupling 
theory |19| ■ This point has been elaborately analyzed in 
an earlier work [I3. In deriving equation [1] we consider 
that the total intermediate scattering function is written 
as, 



0(£)^0a/ctW$>L p W- 



(2) 



Here cf> s hop (t) = s+Kh * {qAF) describes the acti- 
vated hopping dynamics for a single hopping barrier. 
K hop (q,AF) = T{q)P hop {AF). The expression for 
T{q) is given by, T(q) = £(1 - G(q)). For 
the present work the q dependence of T{q) will be 
neglected and F(q) will be set to unity. Ph op {AF) 
is the average hopping rate which is a function of 
the free energy barrier height, AF and is given by 
P hop (AF) = ±exp(-AF/k B T) [H[. The total con- 
tribution from multiple barrier hopping motion is 
written as, ECpW = / <t> s hop {t)V° tatic {&F)dAF = 
f e -tK hop (AF)-pstaticr^p\ ( j i ^p where •p s * a * ic (A_F) 
is considered to be Gaussian, , p s * a * ic (AF) = 

We will call this dis- 



-(AF-AF o )/28AF^ 



V2ir5AF 2 

tribution as the static barrier height distribution. 

We now discuss the MCT part of the calculation. As 
in our earlier model calculation [lj], the values of Q 
and 7 are kept fixed at unity, neglecting their temper- 
ature dependence and the scaling time is taken as lps. 

A i = £ ^+ e (T+(f^ andA 2 = ^ + e (iTTT^W (3- 
The MCT formalism predicts a relationship between A 

and Pmct as Pmct — —log(2)/log(l — X). e is a measure 

of the distance from the ergodic to non-ergodic transition 

temperature of the ideal MCT, T9,, thus e = C T0 . Thus 

c 

to calculate the MCT part of the relaxation we need to 
estimate A and Tj?. These quantities can be calculated 
for systems where the static quantities (like static struc- 
ture factor) are known, but for realistic system due to 
the unavailability of the static quantities the estimation 
becomes difficult. We thus use the following methods to 
estimate A and Tj? . According to the analysis of BR, the 
ergodic to non-ergodic transition takes place at the tem- 
perature, Tj, where the energy landscape properties first 
change [ll|, thus T c ° = T L . According to Sastry et al.^, 
Tl also coincides with the temperature where the stretch- 
ing parameter starts falling. From experimental studies 



we know that for Salol system the stretching parameter 
starts falling at T = 278K (tJ. Thus we estimate that 
T c ° = Tl = 278K. Now for the estimation of A we again 
make use of experimental results. MCT is expected to 
explain the dynamics above T^, thus the stretching pa- 
rameter above Tj? should be equal to 0mct — 0.84 p}. 
We have also mentioned that A and (3mct are related. 
Thus A is fixed in such a way that, above T®, we get the 
correct (3mct- 

Next we discuss the implementation of the hopping 
dynamics for Salol system using the predictions of ran- 
dom first order transition (RFOT) theory. According 
to the RFOT theory, the mean barrier, F(r), which is 
used to calculate the mean barrier height, can be writ- 
ten as F(r) = r r ^ffri[r) " fr 3 Ts c where T K (r) and 
T J 4(r) are the surface energy terms at Tr (Kauzmann 
temperature) and Ta (temperature where hopping bar- 
rier disappears) respectively [17], [2l[ • The temperature 
dependence of the configurational entropy can also be 
given by an empirical formula. [23J s c = s/jt(l — Tk/T) 
where Sfu is a system dependent parameter which is also 
related to the specific heat jump at the Kauzmann tem- 
perature (Acp(T) = Sf lt (T K /T)). For the Salol system 
s flt = 2.65, T K = 175X and T A = 333K. With these 
values of the parameters the mean barrier height and 
the critical nucleus radius are calculated. We find that 
at T = 28QK, the size of the critical nucleus is above 
unity. Thus although Ta — 333-ftT, for all practical pur- 
pose T — 280K should be considered as the tempera- 
ture where activated events start. Note that as found 
in simulations ,15j this temperature is close to Tj?. The 
value of t is fixed in such a way that at T — 280^, 
both the MCT and the hopping dynamics together pre- 
dict a relaxation time which is close to that obtained 
in the experiments [9j. Thus for Salol system we find 
To = 2400ps. The barrier height can also be calculated 
by considering the shape of the nucleating region to be a 
fuzzy sphere [131 ] . The distribution of barrier height arises 
due to the fluctuation in entropy which can be related 
to the specific heat according to the Landau formula, 
< (AS") 2 >= A; B C p ,[24j where AS is the entropy fluctu- 
ation and Cp is the specific heat. This expression can be 
rewritten in terms of configurational entropy fluctuation 
per bead 5s c and heat capacity jump per bead Ac p (T) 

Ss a 



15 



Ss c 



fc B Ac p (T) 



and 



3k B TT K 



-(r*/a) 3 """* <Sc> Y 47T S/lt (T-T K ) 2 0t/a)3- 

Here r* is the droplet radius determined using F(r*) = 0. 
a is the length of the bead. We can also relate the en- 
tropy fluctuation to the width of the barrier height dis- 

tribution, ^E ~ <3 gf c , and thus define the width of 

the static barrier height distribution in terms of entropy 
fluctuation and specific heat. 

Note that from the main equations of motion (eq[T]and 
eqj2]) it is obvious that the hopping dynamics interacts 
with the MCT dynamics and thus effects the MCT dy- 
namical correlation. But, the effect of MCT dynamics on 
the hopping motion is not so obvious. However, we will 
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FIG. 1: Temperature dependence of the time scale of longtime 
a relaxation for Salol. The timescale for the longtime part of 
the total structural relaxation, <f>(t), obtained from experi- 
ments [9( (black circles) and that calculated from the present 
coupled theory, T tota i (red diamonds) are plotted against tem- 
perature. We fit the Ttotai to the MCT phenomenological ex- 
pression, Ttotai ~ [(T - Tl u )/Tl u ] ~ 7 (blue dashed line), we 
obtain T/ 1 ' = 258^. In inset we plot {t MC t - Ttotai) / Ttotai 
as obtained from the theory (red stars). The black dashed 
line is a guide to the eye. The function shows a jump around 
T — 258K. We find that across this temperature the hopping 
dynamics changes its role and below T 1 / 1 * = 258K it plays a 
"direct" role in the dynamics 



show below that the interaction between 4>mct and 4>hop 
generates a renormalized barrier height distribution that 
is explored by the system. A treatment consistent with 
this renormalized barrier height distribution has not been 
attempted, although results should not be too different 
from the ones presented here. 



III. NUMERICAL RESULTS 

We solve equations [1] and O numerically using a nu- 
merical method presented earlier [25| with a minor mod- 
ification [26|. The essential idea involved in the method 
is to separate the slow and the fast variables and treat 
them differently in the convolution. The short time part 
of 4>{t) and 4>Mcr(t) are calculated exactly with very 
small step-size [2T] and they are then used as input to 
carryout the calculation for the long time part of the 
same. The total relaxation time, Ttotai and the stretch- 
ing parameter, fitotai is obtained from the coupled dy- 
namics by fitting the long time part of <p(t) (obtained 
from eqUJ) to a KWW stretched exponential function, 

<t>(t) = Aexp(-{t/ 'Ttotal)f total 

The plot for the relaxation time Ttotai is given in FigJT] 
where we have also shown the experimental results (for 
Salol) and the fit to the MCT phenomenological ex- 



pression. As evident from the figure, the temperature 
dependence of r tota i reveals several interesting physics. 
The Ttotai compares well with the experimental results 
[9| and predicts the correct glass transition temperature, 
T g = 220K. Recall that in our calculation the ergodic to 
non-ergodic transition takes place at T c ° = 278K where 
the activated dynamics also becomes significant. Thus 
we know that the idealized MCT prediction breaks down 
at T = 278K and below this temperature, its validity as 
such is unclear. However when we fit the relaxation time 
to MCT expression (r ~ [(T - r/ i *)/T/ a ] ~ 7 ) we obtain 
a value, T[ %t — 258K, which is in excellent agreement 
with the experimental studies 0, H H3] • 

The advantage of the present scheme of calculation is 
that although it addresses the coupling between the con- 
tinuous diffusion and the hopping motion, one can sep- 
arately analyze their relative contributions to the total 
dynamics and explore the origin of the apparent validity 
of the MCT below T c °. Due to the nonlinear coupling in 
EqsQ] and the activated dynamics plays both a "di- 
rect" and a "hidden" role in the total relaxation. The 
"direct" role is the direct relaxation of <fi via tfihop- How- 
ever, if we analyze the structure of Eq[TJ we find that the 
role of activated dynamics is also to soften the growth of 
the longitudinal viscosity which finally leads to the relax- 
ation of the otherwise frozen 4>mct ■ Thus <fihop plays a 
" hidden" role in the relaxation by helping <Pmct to relax. 
We now analyze these two different roles of the hopping 
dynamics and their effect on the total relaxation. 

Our study shows that both continuous dynamics and 
activated dynamics changes continuously across T[ lt . 
However, when we plot [t M ct - T to tal) / Ttotai against T 
(in the inset of Fig. [TJ we find that this quantity un- 
dergoes a jump like increase in the temperature range of 
the phenomenological MCT transition temperature T/*'. 
This implies that although below T[ %t the continuous dif- 
fusion still remains as the dominant mode of relaxation, 
there is a change of role and increased contribution of the 
activated dynamics. The plot suggest that in the range, 
T[ %t < T < T®, the activated dynamics plays only a 
" hidden" role whereas below T[ %t it also plays a " direct" 
role in the relaxation process. 

Thus in the range, T/ i4 < T < T c °, the role of the acti- 
vated dynamics is primarily to open up the bottle neck to 
start the diffusive dynamics. Once the latter starts, the 
relaxation primarily takes place via the diffusive channel 
and all the MCT predictions hold good. This provides a 
quantitative explanation of the observations reported by 
Kob et al. [12J. The authors find that till close to T/ u 



the dynamics can be described via MCT although needs 
the static inputs to be calculated at a higher effective 
temperature. 

Figure Q] shows the important result that although ide- 
alized MCT breaks down at T c ° = 278^, the functional 
form of MCT continues to describe the dynamical char- 
acteristics till 258K. The latter is a widely known result, 
usually obtained from fittingthe MCT functional forms 
to the experimental data pi, \Q, [jjj] . Thus our study pro- 
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FIG. 2: Temperature dependence of the stretching parame- 
ter for Salol. The stretching parameter for the longtime part 
of the total structural relaxation, <j>(t), obtained from experi- 
ments 6] (black squares) and 7] (black circle) and that calcu- 
lated from the present coupled theory, ptotai, (red diamonds) 
are plotted against temperature. In the same plot we also 
present the stretching parameter predicted by the RFOT the- 
ory using the static barrier height distribution, (blue dashed- 
dot line) and the MCT (green dashed line) (which shows no 
temperature dependence of the stretching parameter) . 



vides an explanation of this intriguing result in terms of 
the interaction between the MCT and activated dynam- 
ics which leads to hopping induced continuous diffusive 
motion. The study further predicts that below T[ %t al- 
though the continuous diffusive dynamics plays an im- 
portant role in the relaxation, the increased contribution 
from the activated dynamics finally leads to the break- 
down of the MCT predictions. 

In Figj2] we plot the calculated temperature depen- 
dence of the stretching parameter 0. In the same figure 
we also plot as predicted by the MCT and the RFOT 
theory (using static barrier height distribution) and also 
the experimental results on Salol [a, 0| ■ It is clear that 
neither the MCT nor the RFOT can alone describe the 
proper temperature dependence of the stretching param- 
eter. However, the coupled theory predicts a tempera- 
ture dependence of fitotai which is qualitatively similar 
to that found in experiments. The stretching param- 
eter is known to provide a measure of the heterogene- 
ity in the system. In our study the heterogeneity arises 
from the distribution of hopping barriers. However, be- 
low T®, f3 'total is smaller than j3rfot which means that 
the barrier height distribution which participates in the 
dynamics (dynamic barrier height distribution) is nar- 
rower than the initial assumption of distribution (static 
barrier height distribution). This modification of the bar- 
rier height distribution (or heterogeneity) takes place due 
to the coupling between the diffusive and activated dy- 
namics. The timescale analysis reveals that the smaller 



barriers contribute to the dynamics, which implies that 
the coupling leads to the facilitation effect. The proper 
temperature dependence of fitotai implies that the theory 
correctly predicts the modified barrier height distribution 
and thus the growth in dynamic heterogeneity. 

To quantify these points, we carry out the following 
analysis. Nonexponential relaxation in the system can be 
considered as arising from superposition of exponentials, 

e -(t/T Kww fKw\v _ 1 e -t/T-pdynamicr r \^ T /g\ 



It is now possible to find the distribution function of the 
relaxation times / p d v namlc (j) by Laplace inverting the 
stretched exponential. Note that it is only possible to find 
■pdynavmc^ w ] ien the individual dynamics is exponen- 
tial. Hence we do this analysis of the distribution of re- 
laxation times (and thus the barrier heights which partic- 
ipate in the dynamics) considering only exponential MCT 
relaxation. Once we obtain the distribution of relax- 
ation times we can convert it into the distribution of bar- 
rier heights using, V dynamic {AF)dAF = V dynamic {T)dT, 
provided r can be expressed in terms of AF. The de- 
tails are presented in the Appendix A. We now analyze 
the dependence oi-p d v namlc (AF) on the MCT relaxation 
time and show that r p d v nam ' lc (AF) not only depends on 
V staUc (AF) but also depends on the MCT relaxation 



time t M ct = K 



-l 

MCT- 



We fix the -p statlc {AF) and change 



the tmct by calculating the MCT dynamics (which is 
coupled to r p statlc (AF)) at different temperatures. The 
results are presented in Figj3] 

At high temperatures, where MCT dynamics is fast, 

V dynarmc( AF ^ j g narrQW (^ j s l arge ) ; anc J t he posi- 
tion of the maximum, AF max , is shifted to smaller AF 
(Figj3]). We also find that as the MCT dynamics slows 
down, the distribution becomes broader and more non- 
Gaussian {(3kww also becomes smaller), larger barriers 
contribute to the dynamics, and AF max shifts to higher- 
values (though it always remains smaller than AF ). 
These findings are in accord with the literature. We know 
that as the temperature is lowered (or the dynamics be- 
comes slower) there is a growth in dynamic heterogeneity. 
Thus we may say that 'P d y namlc i y AF) provides a mean 
field estimation of the dynamic heterogeneity. Note that 
■pdynamic^pj j g a i ways narr0 wer than V static {AF) and 

overlaps with it only in the low barrier side. Xia and 
Wolynes have given a physical interpretation of this ex- 
clusion of the higher barriers from the dynamics, using 
the picture of dynamic mosaic structure Q. Since the ex- 
clusion of higher barriers fastens the hopping dynamics 
thus this is termed as "facilitation effect". Our theo- 
retical model does not explicitly consider these mosaic 
structures. However, we find that the presence of the 
continuous dynamics and its coupling to the activated 
dynamics leads to the facilitation effect. Note that the 
facilitation is strongest when the MCT dynamics is ex- 
ponential and should be weaker when we have stretching 
in the MCT dynamics. 

Next we show that this modification of the barrier 




FIG. 3: Temperature dependence of the probability distri- 
bution of barrier heights which participate in the dynamics, 

pdynamios^py j n thig figure we plot the gtatk barrier height 

distribution, p statlc (AF) which is used in the calculation of 
the total dynamics (black solid curve). From the total dynam- 
ics we again estimate the barrier height distribution which 
actually participates in the dynamics (colored dashed line). 
This has been done for different MCT relaxation time scale by 
chaneine the temDerature of the MCT Dart of the dynamics 
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FIG. 4: Fragility dependence of the stretching parameter. 
The stretching parameter for the total structural relaxation, 
4>(t), obtained from experiments (3 e x P t [J] (black circles) and 
that calculated from the present coupled theory, /3 to tai (red 
solid line) are plotted against the square root of fragility in- 
dex D. As D decreases the fragility of the system increases. 
In the same plot we also present the stretching parameter 
predicted by the full static barrier height distribution (blue 
dashed line) . The experimental and the theoretical values are 
at the glass transition temperature T — T g . 



heights, as predicted by the present theory plays a key 
role in describing the proper relationship between the 
fragility index, D, and the stretching parameter, (3, The 
study of Bohmer et aL[J] predicts a relationship between 
the fragility index D and the _stretching parameter (3 at 
T = T g . As discussed in Ref.Q, Pf^ ic is related to the 
width of the static Gaussian distribution of the barrier 
heights SAF. SAF can be related to the fragility, D , 



SAF 
AF a 



1 . Thus (3hap lc can be related to the fragility 



2VD 



@. Pho a p ic = [1 + {AF /2k B Ty/D) 2 }- 1 / 2 . However, it 
was found that this expression alone (with the full static 
barrier height distribution) does not describe the correct 
relationship between (3 and D Q. 

The details of the calculation are presented in Ap- 
pendix B. In Figj4] wc plot the relationship between 
fitotai and D as predicted by the present theory. In 
the same plot we also show the experimental results 
[J]. We find that the present theory captures the cor- 
rect trend. This is to be contrasted with the predictions 
obtained from the static distribution of barrier heights. 
Note that in all the cases, f3 to tai > Phop lC i w hich implies 
that ■pdynamic^pj is a i ways nar rower than y static (AF). 
Earlier we have shown that r p d y namlc ( K AF) overlaps with 
'P statlc ( K AF) only on the low barrier side. Thus the higher 
barriers do not participate in the dynamics leading to 
facilitation effect. Hence, for a wide range of systems 
the theory predicts a proper modification of the barrier 
heights distribution. Note that it was shown by Xia and 
Wolynes that the static barriers height distribution needs 
a cut-off to describe the relationship between (3 and D [3] • 
The authors gave a physical explanation of this " cutoff' 
using the picture of dynamic mosaic structure. 



IV. CONCLUDING REMARKS 

In this article we present a theory that combines the 
effects of continuous dynamics (through mode coupling 
theory) and the activated dynamics (through random 
first order transition theory) in terms of a unified theory. 
In this unified theory the two dynamics dose not just act 
as parallel channels of relaxation, but they interact with 
each other in a way that in the course of describing the 
dynamics they both modify each others behavior. Cer- 
tain universal aspects of this work arise from this non- 
linear interactions. We show that the unified theory can 
provide satisfactory descri ptio n of the relaxation over the 
whole temperature plane [17} and can explain the tem- 
perature dependence of the relaxation time, r and the 
stretching parameter (3. It also predicts the said relation 
between (3 and D. 

The difficulty of performing MCT calculations for re- 
alistic systems is the non-availability of the static quan- 
tities. This difficulty can be overcome if we can estimate 
the ergodic to non-ergodic transition temperature T®. 
Following the analysis by Burmer and Reichman.[ll| we 
assume that the ergodic to non-ergodic transition takes 



place at the temperature where the landscape energy first 
starts changing, Tl, thus T° = Tl- We next use the ob- 
servations reported by Sastry et al that Tl coincides with 
the emergence of non-exponential relaxation [22| . From 
experimental studies we obtain T? = Tl = 278K [7j. In 
accord with simulation studies, [la] in the present theory, 
the activated dynamics becomes significant, close to T°. 
Thus in our theory, below Tj?, the MCT can only relax 
with the aid of the activated dynamics. In this theory, 
below T° the activated dynamics plays both a "hidden" 
and a "direct" role in the relaxation. The "direct" role 
is when it directly helps in structural relaxation and the 
"hidden" role is when it primarily lowers the longitudi- 
nal viscosity thus helping the otherwise frozen continuous 
diffusive dynamics to be active and help in the structural 
relaxation. The relaxation time T to tai obtained from the 
unified theory agrees well with experiments [9( . Although 
we consider the ergodic to non-ergodic transition and the 
activated dynamics to start at T c °, we find that when 
Ttotai is fitted to MCT expression it shows a MCT tran- 
sition temperature T/ lt = 258K. We find that although 
the continuous and activated dynamics changes continu- 
ously across T/ 1 * the function (tmct -Ttotai) /Ttotai shows 
a jump across T/ 1 '. This implies that hopping dynamics 
changes its role across T/ 1 '. In the range T[ %t < T < T c °, 
it primarily plays a "hidden" role and the dynamics is 
well described by the MCT. This explains the study of 
Kob et al. [12| and also the success of MCT in explain- 
ing experimental results below TJ? [8(. However, below 
T/ J *, there is a sudden increase in the contribution from 
the activated dynamics and it plays a " direct" role in the 
relaxation, which results in the breakdown of the MCT 
predictions . Our study predicts that below T/ 4 ', even 
though the MCT breaks down, the continuous dynamics 
persists and plays an important role in the relaxation. 
Note that in an earlier simulation study we have shown 
the coexistence of the continuous and activated dynam- 
ics and we have further shown that activated dynamics, 
which is usually followed by continuous diffusive dynam- 
ics, relaxes the local stress [14j. Thus the simulations 
studies have already observed the "hidden" role of hop- 
ping which has been predicted in the present study. 

The theory can also predict the correct temperature 
dependence of the stretching parameter. Due to the 
growth of dynamic heterogeneity in the system the (3 
value is expected to decrease with T. While MCT does 
not predict any temperature dependence of the stretching 
parameter and RFOT theory predicts a strong tempera- 
ture dependence of Pf^J 10 , the prediction of the present 
theory is closer to the experimental results. The T de- 
pendence of the Ptotai arises due to the barrier height 
distribution of the activate dynamics. However, we find 
that the barrier height distribution which participates 
in the dynamics r p d v namlc (/S.F) gets modified from the 
static barrier height distribution, r p statlc (AF) which is 
used as an input. The coupling of the MCT dynam- 
ics with the activated dynamics is found to lead to this 
modification. The higher barriers in •p statIC (Ai< 1 ) are 



not present in / p d y namtc (AF), predicting a facilitation 
effect. We have shown that as the relaxation time of 
the system slows down, 'P dynamlc {AF) becomes broader 
and includes higher barriers. Thus we predict that 
■pdynarmc^p^ p r0 vides a mean held estimation of the 
dynamic heterogeneity. We also find that this modifica- 
tion of the barrier height distribution, enables the theory 
to predict a proper relation between the fragility index 
D and the stretching parameter, 0. 
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Appendix A 

According to the present theory, r = -r= tf — 

where K^op = —exp(—AF/kBT). Even for the sim- 
ple exponential MCT relaxation the relationship between 
Kmct and AF is non trivial [26|. However, in the low 
temperature limit we can write Kmct = x a 2 +2\ 'a-i = 
oiKhop, where a — A A -i^ x A _ 1 - A is the Debye- Waller 
factor or the form factor (height of the plateau), thus, 
t = , xttt? — ■ As discussed before, V staUc (AF) is Gaus- 

(a+l)K hop > \ ' 

sian. The mean barrier height, AF , is calculated at 
T = T g , using Salol parameters. The width of the dis- 
tribution, SAF is related to the stretching parameter, 
(^static = [i + (6AF/k B T) 2 }-^ 2 [EIll) and it: is taken 
in such a way that it gives /3£* p* ic = 0.5. With these 
above mentioned parameters we calculate the total re- 
laxation and fit it to a stretched exponential. From the 
fitting we obtain tkww and 0kww ■ This is now fed into 

eqEltO Obtain pdynamic^py 

Appendix B 



Since it is found that the mean barrier height is de- 
pendent on the configurational entropy and almost inde- 
pendent of the system thus for all the systems pjj, we 
calculate AF using Salol parameters (sfu — 2.65, Tk — 
175K, T c ° = 278K). The width of the barrier height 
distribution is fixed according to the fragility index D 
which is taken from Bohmer et al. [4J. This then lets 
us define for each system, the static barrier height distri- 
bution and also the stretching in the hopping dynamics. 
For the Salol system we found that r has to be used 
as a fitting parameter to describe the proper timescale 
of the dynamics. In this calculation we also vary t for 
each system in such a way that at T = T g the relaxation 
time is of the order of 100s. We find that for more fragile 
system we need a larger value of r 

We now need to define the corresponding MCT dy- 
namics for each system. Above the MCT transition tem- 
perature T/ 1 *, the dynamics is expected to be described 
completely by the MCT dynamics. Thus the fitotai above 
T®L should be equal to the 0mct value. For the sys- 
tems like Salol [6| , Glycerol [29| and Silica [3(| where the 
value of the stretching parameters at high temperatures 
are know we use those values to find the MCT param- 
eters. For the systems where the stretching parameter 
at high temperatures are not known we use a reasonable 
value (definitely larger than eX pt{T — T g ) and also vary 
it in small amount along with r to get the proper relax- 
ation timescale at T g . As for the value of Tj? and Tk, we 
keep them same as the Salol system. The exp t values are 
again obtained from Bohmer et al. for the corresponding 
D values which are used to describe the hopping dynam- 
ics. So our only constrain is that the timescale should be 
of the order of 100s at T„. 
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